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We develop a quantitative semiclassical formula for the 
resonant tunneling current through a quantum well in a tilted 
magnetic field. It is shown that the current depends only on 
periodic orbits within the quantum well. The theory explains 
the puzzling evolution of the tunneling spectra near a tilt 
angle of 30° as arising from an exchange bifurcation of the 
relevant periodic orbits. 
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The resonant tunneling diode (RTD) in a magnetic 
field tilted with respect to the tunneling direction, has 
been extensively studied in recent years as a simple exper- 
imental system ^igh manifests the quantum signatures 
of classical chaosErQ. The measured I-V characteristics 
show resonance peaks which evolve in a complex m^ 
ner as magnetic field, 5, and tilt angle, ^, are variec 
The existence and periodicity of these peaks in various 
parameter intervals has been -associated with the exig.- 
tence of certain periodic orbitil and their bifurcationsB'El. 
The link to quantum mechanics has been made by intu- 
itive appeals to Gutzwiller oscillations of the density of 
statesEJ, scaling analyses of the exact quantum spectrumcl, 
and the numerical discovery ofpSequences of wavefunc- 
tions scarred by periodic orbitsB. However, previous to 
this work, it has not been shown that periodic orbits 
indeed determine the quantum tunneling oscillations in 
the semiclassical limit. Indeed, the analogy between the 
tilted well and another well-studied chaotic system, the 
hydrogen slnm in magnetic field, has frequently been 
emphasizedB^; however in diamagnetic hydrogen it has 
been shownB that all orbits closed at the nucleus de- 
termine the absorption spectrum, not just the periodic 
orbits. Moreover, even assuming the importance of pe- 
riodic orbits, the periodic orbit theory of the system is 
sufficiently complex that the contribution of specific peri- 
odic orbits and their bifurcations to the experimental I-V 
data has been controversiaE. Below we derive a quantita- 
tive semiclassical formula for the tunneling current which 
demonstrates that the current is dominated by periodic 
orbits and apply the formula to previously unpublished 
data which reveals an interesting exchange bifurcation in- 
volving four period- two orbits. 

In an RTD under a bias voltage V, tunneling current 
flows from the emitter state through the double-barriers 
confining the quantum well. The data presented are from 
an RTD with a 120nm wide well and experimental details 
are given in Ref. 2. When a large magnetic field (> IT) 



is applied, the emitter state is quantized into the first few 
Landau levels. The electric field is normal to the barriers 
(E = Ez), while the magnetic field is tilted in the y — z 
plane, B = cos^z + sin^y. 

After tunneling into the well through the emitter bar- 
rier, an electron gains large kinetic energy before col- 
liding with the collector barrier. Typically, the electron 
will traverse the well hundreds of times before tunneling 
out and will lose much of its energy due to optic phonon 
emission. Therefore, the tunneling is sequential and the 
resonances are substantially broadened, by h/ropt, where 
the phonon emission time is ~ 0.1 psEJ. For describing 
this limit the jBaxdeen tunneling hamiltonian formalism 
is appropriatetd'tll, which greatly simplifies the problem. 
Using this approach and taking into account that the rel- 
ative barrier widths are chosen so that charge accumula- 
tion in the well is negJLigibleEl, one finds that the current is 
simply proportionall3 to the tunneling rate through the 
emitter barrier, We^w • j = ^e^We^wi where rie is the 
surface concentration in the emitter layer. 

We^w can be calculated from ilip-Fermi Golden Rule 
with the coupling matrix elementE3llil between the wave- 
functions and corresponding to the isolated emii- 
ter and isolated well respectively. In the limit when the 
height of the emitter barrier Uq \^ much larger than the 
injection energy e^, the cyclotron energy TiuOc and the 
voltage drop across the barrier, the coupling matrix ele- 
ment can be simplified to : 



dS^^(x,y,0) 



dz 



(1) 



where the integration is performed over the inner surface 
of the emitter barrier z = 0. 

Due to the translational invariance in the x-direction, 
the classical dynamics within tho-wjell can be reduced 
to two degrees of freedom, 2:,dEI with an effective 
potential V{y,z). The tunneling rate, which is pro- 
portional to the square of the matrix element (|l|), can 
be expressed in terms of the related Green function 
G{yi^zi = 0;?/2,^2 = 0;e). This Green functipn is then 
replaced by its semiclassical approximationliS, which is 
determined by all classical trajectories connecting the 
points (?/i,0); (?/2,0). Defining y = {yi ^y2)/2,Ay = 
yi — y2i it is convenient to introduce the Wigner 
transform of the emitter wavefunction, fw{y^Py) = 
h-^ J dAy^e {y - A^/2, 0) {y + Ay/2, 0) exp {ipyAy/h). 
Since the emitter state involves only the few lowest 
single-particle levels, it can be calculated accurately us- 
ing a variational approach. 
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We then obtain for the oscihatory part of We 

{pi)i {pi)f 

7 

t 



Wosc = j dpy J dy fw {y,Pv) ^^""^(^I^^^ j 



dAy 



X SR. 



exp 



-7 ,_-Si-Py^y 



^opt 



(2) 



where 7 is the index of the classical trajectories, = 
{y — Ay, 0;y -\- Ay, 0; e) is the action integral, ^7r^ 
the traversal time and is the (complex) amplitudell2l. 
The momenta p] and pj correspond to respectively the 
initial and final points of the trajectory 7. The factor 
exp (— t^/ropt) represents the effect of phonon emission 
and suppresses the contributions of trajectories longer 
than Topt (corresponding to ~ 4 bounces with the collec- 
tor). 

The next step depends on the properties of the emit- 
ter wavefunction. If this wavefunction were well-localized 
spatially on the scale of the electron wavelength within 
the well, then effectively it would provide a delta-function 
source and the resulting formula for the tunneling rate 
would involve all closed orbits starting, exactly at the 
injection point, as in the atomic caseB. However, as 
noted above, the emitter state is a linear combination 
of the first few Landau levels and its Wigner function 
fwiViPy) has a typical spatial spread of order the ef- 
fective magnetic length, Ib = \/h/eB cos 0, centered at 
the "injection point" y = yi, and falls off rapidly outside 
this interval. Hence Ay ~ T^/Py ^ Ib ^ V^, and in the 
semiclassical limit ft ^ {) one can expand S^/h to second 
order in powers of Ay near the action of the closed orbit. 
The function Sj{y, y)/h is rapidly varying on this spatial 
scale diS h ^ and, its stationary points correspond to 
the periodic orbitstj. Let y^ be a point of contact with 
the emitter for a given isolated periodic orbit (labelled 
by the index /i) , then all nearby closed orbits can be rep- 
resented to the correct order by quadratic expansion in 
8y = y — y^. li y^ is within Ib (~ fi^^'^) of the injection 
point yi, then this orbit will contribute substantially. In 
this case the quadratic expansion correctly approximates 
the value of {y^y)- li \y^ — yi\ ^ Ib then this expan- 
sion is inaccurate near yi, but the contribution of such 
orbits is negligible due to rapid variation of the phase 
^7 {y^ y)' 0^16 then finds 



Wo. 



16 ^^^^P(-^) 



j dy j dpy 



X fw {y^Py)cos 



^-^^QA^y.Spy) 
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where 



2 {Syf + (m^2 - ^11) ^y^Py - ^12 i^Pyf 



(^11 + ^22 + 2) 

^Py = Py~ {P^Ji)yi the integer is the topological index 
of the periodic orbit, and the 2x2 monodromy matrixE^ 



M = {rriij) is calculated at the emitter barrier. Therefore 
we have shown that the tunneling current depends only 
on the periodic orbits. 

The summation in (^) is performed over all isolated 
periodic orbits, both stable and unstable. Near a stable 
orbit, the classical motion is regular and the trajectories 
are confined to invariant tori in the phase space. One 
can still obtain discrete quantized energy levels from such 
regular motion by insisting that action integrals over the 
torus be quantized. This leads to a sequence of ei^ 
functions localized on the tori of the stable islandstSI 
In contrast, the classical motion near an unstable orbit 
is chaotic with no locally conserved actions. As is now 
well-known, such motion cannot be semiclassically quan- 
tized to yield discrete energy levels. Therefore one would 
expect a qualitative difference between the contributions 
of stable and unstable orbits in Eq. (||). This difference, 
can be displayed explicitly by performing the summation 
over repetitions of the primitive periodic orbits. This 
summation can be performed exactly, yielding: 



X I dy I dpyf^{y,py)g^'^ {y,py) 



(4) 



where A (cr, p) 



sinh(cr) 



and the index + or 



cosh(cr)— cos(p) ' 

denotes stable or unstable orbits. The quantity ^/r^ 
is an effective level-broadening which differs in the two 
cases. 

(i) Stable orbits Here we expect the level-broadening 
to arise only due to phonon emission and indeed we find 
r^^ = Topt- If phonon scattering is absent, 00. 
Noting that limcr^o A (a, p) = 27r5 (p), one finds that the 
stable orbit contribution to Wosc is proportional to the 
sum S{e—en,£)' The energies En/ are determined from 
the semiclassical quantization condition of the longitudi- 
nal energy = e — huj^~^{£ + 1/2) along the periodic 
orbit Sjj,{sn,i) = 27rfi(n+^). Due to the harmonic 
approximation the quantization of the transverse oscilla- 
tions around the PO simply yields equally spaced levels 
with spacing huj^~^ , where th^ frequency uj^~^ = 
with (/)^ the winding numbertj and the period of the 
orbit. We find that the coefficient functions in Eq. are 
the Wigner transforms of the harmonic oscillator wave- 
functions corresponding to these transverse modes: 

g';'^{y,Py) =(-l)%(2|QJ)exp(-|QJ) (5) 

where L£ is Laguerre polynomial and is defined in (||). 
Since the result (||) is based on the harmonic approxima- 
tion within a stable island, it is only valid for quantum 
numbers i such that the oscillations remain within the 
island. Equivalently, the number of such modes ^max is 
given by the ratio of the island area to h. The effect of 
phonon scattering according to (H) is to smear out each 
of these delta-function contributions to Wosc over a scale 

^/Topt- 
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(ii) Unstable orbits. For unstable periodic orbits we 
find that the "effective" relaxation rate 1 / r^^ = 1/ Topt + 
(-£ + 1/2) A^, where is the Lyapunov exponent near the 
orbit /i. Hence Teff is finite and equal to the Lyapunov 
time when Topt oo. Therefore, instability acts as a 
sort of intrinsic level-broadening and the contribution of 
unstable POs to never approaches a delta function 
describing individual levels. Instead this contribution de- 
scribes the well-known clustering of levels responsible for 
Gutzwiller oscillations of the density of statesEj, and in 
the tunneling rate each PO contribution is weighted by 
the function 



.sin {S^/h - 7rn^/2) 



1 



)} 



(6) 



sinh (T^/<ff) 

corresponding to Wigner functions average(E^ over the 
eigenstates of the cluster. 

It is not clear in this case how to introduce the cut- 
off value ^max, however for each PO the contributions to 
Wosc decay exponentially ~ exp {—X^T^i) with increas- 
ing and we find that for our case the main contribution 
to the tunneling rate is given by the £ = term. 

Above we have argued that among the POs which 
reach the emitter only those with points of contact 
within a distance ~ Ib from the injection point con- 
tribute strongly. Our formulas (4), (5), (6) permit us 
now to propose a precise criterion : a PO is "accessi- 
ble" to tunneling electrons if the emitter Wigner func- 
tion centered on y = py = and of spatial width 
~ Ib overlaps the functions g^'^ centered at the contact 
points y^ji^{p^)y' The localization length of the wave- 
functions corresponding to the functions can 
be sho wn by an extension of our above analysis to be 

= ^J2n\m}[^\l ^\^-Tt^ [M^]|. The spatial scale as- 
sociated with an unstable orbit is of the same form but 
instead of Gaussian decay of the wavefunctions at ^y > 
one finds rapid oscillations ~ exp (i(5?/^/2/^) . 

The relations (§) , (|5|) , (^) constitute a precise semiclas- 
sical expression for the tunneling rate in terms of the 
contributions of the distinct periodic orbits which reach 
the emitter barrier ("emitter orbits"). We will apply 
them now to a specific parameter regime in the tilted 
well to compare qualitatively and quantitatively to the 
experimentally-observed I-V characteristics. 

In the recent periodic orbit theoryS is was shown that 
within the set of periodic orbits (POs) which collide with 
the collector n times (period-n orbits), there exist orbits 
which collide with the emitter m times, where m < 
and it is useful to classify POs by the two integers (m, n). 
At ^ = the only resonances observed in the I-V charac- 
teristic are associated with Bohr-Sommerfeld quantiza- 
tion of the (1,1) orbit which traverses the well with zero 
cyclotron energy. When such period-one "travers- 
ing" orbits still exist in some regions of the B — V pa- 
rameter space and define a background frequency of reso- 
nance peaks. However now additional resonances appear 



corresponding to doubling or tripling of the frequency 
of peakad. These new peaks are associated with the exis- 
tence of period-two and period-three orbits which appear 
and disappear as a result of bifurcationsO. Here we fo- 
cus on the peak-doubling in the interval 29° < 6 < 34°, 
where, as was shown in Ref. 5, there are four most rel- 
evant orbits, denoted by (l,2)i, (1,2)2, (1,2)^, (1,2)^. 
For 6 < 30°, the starred and unstarred orbits are paired 
in the sense that they eventually become identical and 
disappear in inverse tangent bifurcations. 

The evolution of these orbits is represented by the four 
colored lines in Fig. 1. We recallB~0 that under exper- 
imental conditions the classical mechanics depends only 
on two parameters /3 = 2Bvq/E (where vq is the veloc- 
ity corresponding to the total energy : vq = ^y2e/m*) 
and on the tilt angle 0. As /3 increases from zjejfou-these 
four (1,2) orbits appear in cusp bifurcationsOot^ and 
then disappear pairwise at higher p in the inverse tan- 
gent bifurcations already mentioned. The hatched region 
denotes the semiclassical width of the emitter state, while 
the gray-scale regions denote the semiclassical widths de- 
fined by the (evaluated here at 5 = 8 T) for the 
most accessible of these orbits for a given value of ^tZl. 
Wherever these regions overlap the semiclassical formula 
(4) will predict peak-doubling regions to appear in the 
B — V parameter space for the corresponding values of 

as seen in Fig. 2. 

A fascinating feature of the classical dynamics, noted 
in Ref. 5, occurs near = 30° for the experimental pa- 
rameters of Ref. 2.|-The four (1,2) orbits undergo an 
exchange bifurcationll3 so that for > 30° the (l,2)i is 
paired with (1, 2)2, whereas the (1,2)^ orbit is now paired 
with the (1,2)2. For ^ 29° the (1,2)2 orbit is accessi- 
ble for over its entire interval of existence, 4.3 < P < 8. 
and the (1,2)^ is accessible in the overlapping interval 
7.9 < /3 < 10.9. This situation manifests itself in a large 
and continuous region of peak-doubling in the experi- 
mental dataH. By ^ = 31°, the pairing of the (1, 2) orbits 
has been interchanged (Fig. lb); the (1, 2)^ — (1, 2)2 pair 
exists at lower /3 (higher voltage), and is moving away 
from the semiclassical accessibility region, whereas the 
(1,2)1 orbit has become most accessible. Thus we ex- 
pect the peak-doubling region to split into two smaller 
regions (Fig. 2b) and the amplitude of the oscillations 
to be weaker in the low-magnetic-field, high-voltage re- 
gion (Fig. 2a). By 6 = 34°, the (1,2)2 orbit has become 
inaccessible and these oscillations are no longer seen in 
Figs. 2d and 2e. In contrast, the (1, 2)i orbit remains ac- 
cessible and still produces strong oscillations in the data 
(Figs. 2d,e). 

In Figs. 2c and 2f, we show the results of the semi- 
classical theory at = 31° and = 34°. The semiclas- 
sical calculation of the tunneling current using formulae 
is in good agreement with the experimental 
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FIG. 1. (color) Bifurcation diagrams for the relevant period-two orbits at (a) before (Q — 29°) and (b) after (Q — 31° and 
34°) the exchange bifurcation. The shading represents the calculated localization lengths of the wavefunctions localized near 
the most relevant (1,2) periodic orbits (see text). The hatched region denotes the semiclassical width of the emitter state. 
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FIG. 2. (a),(d) Samples of experimental resonant tunneling 1 — V traces, from which peak positions are determined and 
(b),(e) plotted versus total magnetic field, to be compared with (c),(f) results of the semiclassical calculation. Note at (a-c) 
^ = 31° the onset of the separation of the large region of peak-doubling immediately after the exchange bifurcation and, at 
(d-f) Q — 34°, the peak-doubled region at lower magnetic field created by the exchange bifurcation is almost invisible due to 
much weaker coupling of the corresponding orbit (1,2)2 to the emitter (see Fig. Ic). The I-V traces (a), (d) evidence this by 
the disappearance of the oscillations at higher bias voltages. 



5 



